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Quantum critical (QC) phase transitions generally lead to the absence of quasiparticles. The 
resulting correlated quantum fluid, when thermally excited, displays rich universal dynamics. We 
establish non-perturbative constraints on the linear-response dynamics of conformal QC systems 
at finite temperature, in spatial dimensions above one. Specifically, we analyze the large fre¬ 
quency/momentum asymptotics of observables, which we use to derive powerful sum rules and 
inequalities. The general results are applied to the 0{N) Wilson-Fisher fixed point, describing the 
QC Ising model when = 1. We focus on the order parameter and scalar susceptibilities, and the 
dynamical shear viscosity. Connections to simulations, experiments and gauge theories are made. 


The quantum Ising model in two spatial dimensions 
(2-l-lD), e.g. on a square lattice, undergoes a quantum 
critical (QC) phase transition as the ratio of the trans¬ 
verse magnetic field to the exchange coupling is tuned. 
It is the archetypal example of a non-trivial 2-l-lD QC 
point, possibly the simplest one with Z 2 symmetry, but 
lacks an exact solution contrary to its lower dimen¬ 
sional counterpart. Rather than having quasiparticles 
excitations, present in the para/ferromagnetic phases, 
the spectrum at the QC point is continuous. Various 
methods such as Monte Carlo simulations,^ field theory 
expansions,and recently conformal bootstrap,^ have 
shed light on the critical exponents characterizing its 
thermodynamics and groundstate correlations. In con¬ 
trast, little is known about its quantum dynamical prop- 
ertiesx at finite temperature,^’^ which are not only impor¬ 
tant to understand the nature of this strongly correlated 
quantum fluid but also of clear relevance to experiments. 

In this article we study QC dynamics, with a focus on 
the quantum 0{N) Wilson-Fisher fixed point which de¬ 
scribes the QC transition for the quantum Ising {N = 1) 
and XY {N = 2) models, and the Neel transition in cer¬ 
tain antiferromagnets {N = 3). Focusing on a large class 
of experimentally relevant observables, we establish non- 
perturbative results for the large frequency/momentum 
asymptotic behavior and sum rules. These provide strong 
constraints on the universal scaling functions character¬ 
izing the system’s low-energy responses. The exact sum 
rules can be seen as generalizations of the celebrated /- 
sum rule to scale invariant systems. Our results provide 
rigorous means to assess approximations, constrain nu¬ 
merical results, and ultimately assist with the analysis of 
experimental data. The methods we use partly rely on 
the conformal symmetry of the QC point, present for the 
0{N) Wilson-Fisher fixed point. However, the key ideas 
are more general, and they greatly generalize the recent 
analysis® for the dynamical conductivity of 2-|-ID confor¬ 
mal field theories (CFTs). The paper is organized as fol¬ 
lows: We first establish general properties regarding the 
asymptotics and sum rules of CFTs, and subsequently 
apply them to the Wilson-Fisher theory, and finally give 
a broad outlook, including a discussion regarding the im¬ 
plications for Monte Carlo simulations. 
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FIG. 1. a) Phase diagram near a quantum critical point 
(QCP). b) Asymptotic behavior of the Euclidean suscepti¬ 
bility associated with an operator O of scaling dimension A: 
X^iioJn) = {0{—uJn)0{ujn))T. c) Schematic operator product 
expansion (OPE) determining the asymptotics of y. “desc.” 
denotes the descendants of the primary On (dimension A„). 


Asymptotics and OPE: We consider a thermally 
excited system tuned to a QC point via a non-thermal 
parameter g. In the phase diagram Fig. la, this cor¬ 
responds to the line in the QC fan at g = gc and 
T > 0. We are interested in the linear-response dy¬ 
namics at finite temperature, more precisely in the re¬ 
tarded dynamical susceptibility associated with a bosonic 
observable O, such as the energy or charge density: 
X^(t,x) = —z0(t) ([C>(t, x), 0(0,0)]).p, where the aver¬ 
age is taken over the thermal ensemble. We set h = 
fcs = c = 1; c is the characteristic speed near the QC 
point. We will often work in Fourier space: x^(w,k) = 
/dM‘^xx'^(t, where w is the real frequency 

and k the momentum. Using T, the only energy scale 
available, can be rewritten to make its scaling prop¬ 
erties manifest: 

, ( 1 ) 

where Aq is the scaling dimension of O, d) is a universal 
scaling (complex) function, and z the dynamical critical 
exponent. This scaling structure emerges at low energies. 
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i.e. u!, |k|,r <C Auvi where Ayv is a microscopic lattice 
energy scale, represented by the horizontal dot-dashed 
line in Fig. la. We emphasize that in this regime the 
ratios (w, |k|^)/T can be arbitrary. We introduce the 
corresponding universal response function 


TZ{oj, k) 


X^(^,k) 

ioj - 0+ ’ 


( 2 ) 


using the Kubo prescription. E.g. if O = J^, is a con¬ 
served current, Xxx i® tii® a:a;-polarization function and 
7Zxx(!^,0) = axx(oj) the dynamical conductivity. Due to 
the strong interactions and the resulting absence of quasi¬ 
particles in generic QC systems, little is known about 
these universal responses, and our goal is to unravel some 
of their robust properties. First, let us begin with the 
large-frequency regime, uj ^ T, |k|, where the dynamics 
are near those of the groundstate. These can be elegantly 
studied via the operator product expansion^ (OPE) of O 
with itself. The OPE is an operator relation and does 
not depend on temperature. For a general QFT, it is a 
short time/distance expansion that captures the behav¬ 
ior of the operator product Oi{t, x)O 2 ( 0 , 0 ) as t, |x| —>• 0, 
which by locality, can be expressed as an infinite sum of 
operators evaluated at t, |x| = 0. We will mostly focus 
on CFTs, which have z = 1 and describe a large class 
of experimentally relevant QC phase transitions such as 
those in the quantum Ising and XY models. In a CFT, 
the OO OPE®’9 of a primary operator O with scaling 
dimension Aq reads (Fig. Ic) 
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stituting this into Eq. (4) we obtain a general expression 
for the large-fc asymptotic expansion of y. We see that 
the lowest dimension operators appearing in the OPE 
dictate how the susceptibility approaches its groundstate 
value as T/|fc| —>■ 0. To obtain the real quantum dy¬ 
namics, we can analytically continue the imaginary fre¬ 
quency expansion Eq. (4) to real frequencies^*^ termwise, 
with the replacement iojn oj + iO’*". This follows from 
the structure of the OPE and the spectral representa¬ 
tion connecting the Euclidean and retarded susceptibili¬ 
ties (see App. A for an extension of the proof in [11]). 

Interestingly, unitarity and conformal symmetry con¬ 
strain the scaling dimensions of these operators:A„ > 
{D — 2)/2. This leads to important inequalities for dy¬ 
namical susceptibilities. Let us work in 2-l-lD and con¬ 
sider a putative low-energy susceptibility x(w) obtained 
from an experiment or simulation, and express it as 


X(w) 


, ,2(5i-3 
— UJ 



(5) 


at large frequencies. Finding either 5i < 1/2 would vio¬ 
late unitarity bounds and thus rule out a conformal QC 
point.For Wilson-Fisher QC points, we shall see that 
the stronger condition, ^2 > 1.4, holds. Before applying 
the above general results to those CFTs, we discuss how 
the asymptotics can be used to prove sum rules for any 
susceptibility (2-point function). 

Sum rules: We put forth a powerful sum rule for the 
real frequency quantum dynamical response function: 

/ OO 7 

— Re(57^(a;,k) = —. (6) 

-OO 'X 


which is expressed in imaginary time r: |a;p = r^-|-x^. A 
primary operator transforms homogeneously under con¬ 
formal transformations; e.g. conserved currents and the 
order parameter in the 0(A) model. The sum in Eq. (3) 
is over primaries On with scaling dimensions A„; it in¬ 
cludes the identity (dimension 0). The differential opera¬ 
tor Cn{x, is homogeneous under x — bx, and encodes 
the contributions from the descendants of On (obtained 
by applying derivatives to 0„). Going to Fourier space 
and taking a thermal expectation value (TEV) we ob¬ 
tain a key result (Fig. lb): the jfcj ^ T behavior of the 
Euclidean susceptibility, 

yE(fc) = ^ E(A)i^ + ...V (4) 

where jfcp = -|- k^, = 27rTn is a Matsubara fre¬ 

quency. The dimensionless functions c„ encode the ap¬ 
propriate fc-space tensor structure (and can contain log¬ 
arithms). The dots correspond to higher powers of r/jfcj 
arising from the descendants of On. Crucially, a scaling 
operator will acquire a TEV, {On)T = dn , since T is 
the only energy scale, dn is a universal real number. Sub- 


5TZ is defined as in Eq. (2), with a modified susceptibil¬ 
ity defined below. The sum rule is independent 

of small frequency details, and fundamentally relies on 
the retarded causal structure of y^. More precisely it 
is the zero-frequency limit of the Kramers-Kronig trans¬ 
form for the modified susceptibility, 5y^(tLi,k) — ^y°°, 
which we now discuss. In the QC scaling regime, y 
does not usually decay at large frequencies unlike on the 
lattice because it encodes excitations at all scales. To 
formulate the sum rule, we thus generally need to sub¬ 
tract terms, denoted by y, from y to remove its large-w 
divergence: ^'*“**^ ^y(D,k) = y(r2,k) — y(D), where D is 
a complex frequency in the upper half-plane. In some 
cases, one further needs to subtract a remaining constant: 
(5y°° = ^y^(D—>- 100 , k), where the limit is taken at fixed 
T. We emphasize that y(D) is momentum-independent 
because the asymptotic behavior Eq. (4) depends on pow¬ 
ers of Tjyjojn + jkp due to the asymptotic re-emergence 
of Lorentz invariance (broken by T), and we fix k as we 
take uJn 3> T. We note that the correlation functions 
studied here only depend on the magnitude of the mo¬ 
mentum jkj, which implies that Imy^,Re7?. are w-even 
functions, so that the integral Eq. (6) can be written for 
w > 0. 
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The highly non-trivial and theory dependent informa¬ 
tion is contained in the subtraction terms y, which 
are determined from the large-frequency behavior, i.e. 
from the leading terms in OPE, Eq. (3). We now derive 
some general properties of the subtractions. First, the 
main subtraction (5% = x ~ X is generally required be¬ 
cause the leading asymptotic |fc| T behavior of x^(^) 
is , and most operators have 2IS.o > D. In con¬ 

trast, in almost all cases the subtraction of a constant 
is not needed, i.e. 6x°° = 0. Indeed, from Eq. (4) this 
constant can be non-zero only if the OO OPE contains 
an operator O* with dimension A* = 2Ao — D. More¬ 
over, Ot needs to have a non-zero TEV. In which case, 
5x°° oc {OAjT and the constant of proportionality is the 
corresponding OPE coefficient. A further necessary con¬ 
dition for 5x°° 7 ^ 0 is Ao > (3II — 2)/4 because of the 
unitarity bound^^ on A*. A generic case where the sub¬ 
traction Sx°° appears is for a 2-point function of the 
stress tensor, because the latter has scaling dimension D. 
In this case O* = since the stress tensor generally ap¬ 
pears in the OPE. Below we will the consequences 

of this for the shear viscosity. 

0{N) model: We now apply the above general re¬ 
sults to the QC point of the quantum 0(A) modeP’^® 
in dimensions 2 < D < 4. This is the famous Wilson- 
Fisher conformal fixed point. It describes a variety of ex¬ 
perimentally relevant quantum phase transitions: Ising 
{N = 1), XY (A = 2), etc. An exact solution exists at 
A = oo, which we will use to perform non-trivial checks. 
As a field theory, the 0(A) (non-linear sigma) model 
is defined by the action S = Jd^x^df^(padfj,(pa, where 
ipaix) is a real A-component vector field of fixed norm 
V’a^fia = 1. As the coupling g is increased the system un¬ 
dergoes a QC phase transition at g = gc from a broken 
symmetry phase to a symmetric one for g > gc (Fig. la). 

For our asymptotics/sum rule analysis we need 
the list of operators (On, An) with low dimensions 
An < D. These are known from large-A and small 
(4 — D) expansions,Monte Carlo, ^ non-perturbative 
bootstrap,etc. The first one being the order pa¬ 
rameter field (j)a with dimension A^ = (D — 2 + rj^)/2, 
where rj^ is the held’s anomalous dimension. The follow¬ 
ing 0(A)-invariant operators will also appear: the “ther¬ 
mal” operator (Og, Ag), the conserved currents — 

1), and the stress tensor (Tgn,D). The dimensions of 
the currents and stress tensor receive no anomalous cor¬ 
rections because they are protected by symmetries. The 
operator Og (often denoted by e in the context of the 
Ising model) is associated with the Lagrange multiplier 
field A(a;) that constrains (pa^a = 1 in the 0(A) model. 
It has dimension Ag = D—ljv, where v is the correlation 
length exponent; for the D = 3 Ising case,^’"^ Ag = 1.413. 
It is directly related to the singlet (jp, and tunes the sys¬ 
tem away from the QC point. Being the only relevant 
0(A)-symmetric scalar, it is the most important opera¬ 
tor as it dominates the asymptotic quantum dynamics: 
we will see that it generally gives the first finite-T correc¬ 
tion. This was recently shown® to be the case for the con¬ 


ductivity of the 0(A) model, and observed numerically® 
for A = 2. Given the generality of our OPE analysis, we 
infer that this “dominance” of the relevant symmetric 
scalar is a generic property of QC transitions. 

Order parameter susceptibility: We hrst study 
Xab{k) = {4>a(—k)(l)b(k))T, i.c. the order parameter sus¬ 
ceptibility. It is one of the simplest observables, and 
yields the low-energy staggered spin susceptibility of 
quantum antiferromagnets with transitions in the 0(A) 
universality class. We begin by analyzing its asymptotics. 
By symmetry, and from the knowledge of the operators 
with low dimensions we can write the leading terms in 
the (jiatl^b OPE: 






Oct><j>gOg(d) 

^2A^-Ag 


Ctptpxx fj,Xi/Tg^j^(0) 

^2A^ —Z)+2 


(7) 


where we focus on a, & = 1 since Xab is diagonal by virtue 
of 0(A) symmetry. We have omitted the contribution 
from the currents because they have vanishing TEV 
(no excess charge or net current in the thermal ensemble). 
Taking the TEV of Eq. (7) gives the asymptotic behavior 


2A*-D 


Xii(ia;„,k) = |A:| * 
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C^/C^ in Eqs. (7)/(8) are real OPE coefficients in po¬ 
sition/momentum space, which can be obtained from 
groundstate 3-point functions. As anticipated, the first 
subleading term comes from the relevant scalar Og. The 
next term arises from the stress tensor, where {Tgn)T = 
is diagonal. At A = OO, A^ = {D — 2)j2 saturates 
the unitarity bound, but finite A fluctuations lead to a 
small anomalous dimension <C l.^’^dds OPE co¬ 
efficients C^^g , Ccfnf,T are generally finite and can be com¬ 
puted using a 1/A expansion for instance. C^^g has 
been computed using bootstrap^^ and Monte Carlo^^ for 
A= 1. From the above expansion, we can derive the sum 
rule for Xab- First, for any A, Xab decays sufficiently 
fast at large frequencies so that the subtractions vanish, 
Xab = dx^b — takes its simplest form: 


poo 

/ dwRe7^a&(w, k) = 0, 
Jo 


(9) 


where 77.ab(w,k) = Xa&(‘^)l^)/(*^ ~ is the response. 

When A = OO, we have the exact solution for 2 < 
D < 4: Xab(*Wn,k) = Sabliujl -I- k^ -I- mf.), where 
iriT = ©dT is the thermal mass, and 0^ is a positive 
number,^® see App.B. Expanding for |fc| T, we get 
Xab(fc) = ^ [1 - + (^)'‘ + •■•]■ In agreement 

with the OPE, Eq. (7), the subleading term —m'^/k'^ has 
A^^°° — 2 (i.e. 1/v = D — 2) and is proportional to 

{Og)T = \/Nrn\. This later TEV is evaluated® in the 
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N = oo limit. We note the absence of a contribution 
from the stress tensor, ^ m^/\k\^. Although the real- 
space OPE coefficient in Eq- (J) is non-zero, upon 

Fourier transforming to A:-space, that term does not con¬ 
tribute to the \arge-k behavior. This is an artefact of 
N = oo, where (j) has no anomalous dimension. Finally, 
the sum rule Eq. (9) can be easily checked as the spectral 
function is a sum of (quasiparticle) delta functions. 

Scalar susceptibility: The scalar susceptibility Xs 
is the 2-point function of the “thermal” operator, 
{Og{—k)Og{k))T- It has recently been the focus of atten¬ 
tion in the study of the amplitude “Higgs” mode.^"^^^® 
Again, we first examine the OgOg OPE. The terms rele¬ 
vant here are given, mutatis mutandis, by Eq. (7). This 
then leads to the large-fc expansion Eq. (8) with (0^, A^) 
replaced by {Og,Ag). With this data, we can derive the 
sum rule for Xs- First, Sx"^ = 0 since there is no 0(A)- 
singlet with dimension A* = 2Ag — D in the spectrum. 
The other ingredient needed to build the sum rule is 
the term removing the large-w divergence, Xs- In this 
case, it is simply the groundstate value of Xs at k = 0: 
Xs(fl) = 0) = reads: 

dwRe[7^,(a;,k)-7^f=°(a;,0)] = 0. (10) 

We can again carry out the asymptotic analysis ex¬ 
actly for N ^ oo. The result is (App. B): xf (fc) = 

- Sif P - -)> ■^"here are 

A-dependent constants. Interestingly, the coefficient of 
the subleading term, Og, vanishes exactly for D = 3. This 
comes from the somewhat surprising fact that the ther¬ 
mal operator Og does not appear by itself in the OgOg 
OPE when I? = 3 in the N = oo limit. In other words, 
the OPE coefficient vanishes. This does not hap¬ 

pen for D ^ 3, and we do not expect it to hold at finite 
N in D = 3. Indeed, for the Ising case this coefficient 
was recently computed using Monte Carlo methods and 
found to be finite.®® Finally, the sum rule Eq. (10) can be 
checked numerically at A = oo (App.B). 

Dynamical shear viscosity: Einally we examine a 
correlator involving the stress tensor. Not only is this 
of fundamental interest because it can be defined for 
any GET, but it will also reveal the full complexity of 
the sum rule. We consider the dynamical shear viscos¬ 
ity, r]{u},k) = Xg (^^,k)/(*a; — 0+), obtained from the 
Txy 2-point function, Xy- 'I'xy measures the flux of x- 
momentum in the y-direction, and t] probes the system’s 
resistance against momentum gradients. The asymptotic 
behavior of r] follows from the T^yT^y OPE, which we 
here formulate in momentum space: 

lim T^y{k)T,cy{-k + p) = CT\k\°S{p) 

|fc|>|p| 

+ CTTg\k\^-^^Og{p) + + ■■■, ( 11 ) 

where here k = 0 for simplicity. This can then be used 
to derive a sum rule for r], which is more involved than 



for the response functions considered above. Eor one, 
(5x“ = is non-zero, as was explained above 

on general grounds for 2-point functions involving 
Second, the subtraction involved in Sxy is temperature 
dependent because Og is relevant. This leads to the fol¬ 
lowing sum rule for the shear response: 

pOO 

/ dujlie[r]{uj,'k) — = CyP, (12) 

Jo 

where d = D - 1, A = CTTg{Og)T, P = {Txx)t is the 
pressure of the GET, and Cg = —ttSx'J^ /(2P) is a dimen¬ 
sionless constant. The second term in the integrand is 
, 0 ) and mirrors the subtraction in the scalar sum 
rule. The third one depends on temperature via A oc 
and scales with a non-trivial u power depending on the 
correlation length exponent v via Ag = D — Xjv. Some 
QG theories are simpler in that they lack a relevant scalar 
that condenses at T > 0, as we now discuss. 

We contrast the above shear sum rule with the simpler 
ones obtained^^’^® for M = d super Yang-Mills and pure 
Yang-Mills, which are gauge theories in D = A. In those 
cases, the result is as in Eq. (12) except that the third 
term in the integrand is absent. This stems from the fact 
that those theories do not contain a symmetric relevant 
scalar like Og, i.e. they are not realized by fine tuning a 
symmetric “mass” term. The massless version of QED in 
D = 3 with many Dirac fermions coupled to a U(l) gauge 
field also satisfies this property, being a stable phase. It 
will thus have a shear sum rule of the same form as super 
Yang-Mills. Finally, we note that shear sum rules anal¬ 
ogous to Eq. (12) were derived in the context of strongly 
interacting ultracold Fermi gases,which generally do 
not have emergent Lorentz symmetry. 

Outlook: Our non-perturbative results, via the op¬ 
erator product expansion (OPE), for the asymptotics 
and sum rules apply to a wide class of conformal QG 
points, many of which describe experimentally relevant 
systems. It will be interesting to apply the program de¬ 
scribed in this article to theories other than to the 0(A) 
Wilson-Fisher fixed point, treated here, or even to non- 
conformal QG systems. The strong constraints we have 
derived will also be useful for the analysis of numerical 
and experimental data. Eor instance, quantum Monte 
Garlo is a powerful tool to study QG dynamics in imagi¬ 
nary time,®’®^^®* and can be used to study the asymptotic 
regime where the OPE analysis applies, as was recently 
shown® for the conductivity. The asymptotics and sum 
rules will also help with the difficult task of analytically 
continuing the imaginary time data to real time by con¬ 
straining the allowed scaling functions. Along those lines, 
our results can be used with a novel method®’®® of an¬ 
alytic continuation based on the AdS/GFT holographic 
principle:®® Specific data about a QG theory can be en¬ 
coded in holographic physically-motivated Ansatzes for 
the scaling functions. These can then be used to perform 
the continuation. 
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Appendix A: Analytic continuation of asymptotic susceptibility 

This appendix explains how the quantum critical (QG) dynamics in the near groundstate regime follow from the 
operator product expansion (OPE) in imaginary time. More precisely, we show that the asymptotics at large real 
frequencies, w ^ T, can be obtained from the Euclidean frequency result via term by term analytic continuation. To 
do so, we adapt the line of reasoning put forth in Ref. 11. The analysis begins with the spectral representation for the 
Euclidean susceptibility: 



dw' p{uj') 

TT uj' — iU,' 


(Al) 


where p{uj) = Imy^(a;+j0^) is the spectral density of the retarded 2-point function. We have dropped the momentum 
dependence of p and any other indices as these are not crucial to the discussion. Here 11 > 0 so that Al is a positive 
frequency along the imaginary axis, so that taking a principal value of the integral is not necessary. At finite 
temperature we should strictly speaking use Matsubara frequencies, H = a;„ = 27rrn, in Eq. (Al). However, it is more 
convenient to employ the spectral representation to analytically continue the susceptibility to arbitrary frequencies 
along the imaginary axis, and more generally to the upper half of the complex plane. In writing the above equation, 
we have assumed that p{uj) vanishes as |a;| —?► oo. If this fails one needs to use a subtracted density, 5p = Im^y^, 
as was discussed in the main text. Our argument is independent of this complication, on which we shall comment 
towards the end of the section. 

For the cases of interest, namely 2-point functions of bosonic operators, p{ui) is w-odd allowing us to reduce the 
integral to positive frequencies: 



duj' 2uj' 

TT UJ'^ + 


P(w') • 


(A2) 


We use the convention according to which p{uj) is positive for w > 0. Now, assuming that p{uj) has a large-frequency 
expansion in powers of I/w, we will show that it leads to a corresponding expansion in powers of I/H for y^(if2). 
Logarithms can appear but do not spoil the correspondence, see Sec. Al. Let us assume that p{uj) is bounded by 
the tail l/w'^, 0 < e < 1, at large frequencies. Using the positivity of Eq. (A2), one can show that is bounded 
by l/H*^ as H —>■ oo, and thus vanishes in that limit. It can also be shown that the first derivative of Eq. (A2) with 
respect to 1/H diverges as H —>■ oo. This establishes that y® —>■ l/fl*^ at large H, with an exponent 0 < e < e. The 
procedure can be adapted to a more general spectral tail, where £ is a positive integer. In this case, one needs 

to take £ -I- 1 derivatives with respect to I/H. The argument can be further iterated for all the terms in the power 
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law expansion of p by considering a modified density p from which one has subtracted terms up to the one that is 
targeted. We have thus established that an expansion for p in powers of l/w leads to a corresponding expansion in 
powers for y®. We now turn to the crux of the proof and show that these expansions are precisely related by termwise 
analytic continuation, e-)- a; + iO’*'. 

Let us assume that p{uj) contains the power law l/w“ for cu ^ cuo, where we are mainly interested in the case where 
the infrared scale is temperature, ojq = T. Consider a > 0 but not an even integer, then using Eq. (A2) we get 

/■“ duj' 2J 1 _ 1 2 ^ 

TT uj'°‘ sin(7ra/2) ^ tt ’ 

where a„ = (—l)"+^/(a — 2n). We see that the first term is the expected power law, and is independent of the 
infrared scale ujq defining the asymptotic regime. The remaining terms depend on wg and form a series that contains 
only even powers of l/fl. Due to their incompatible powers, these do not interfere with the first term. In addition, 
the factor of l/sin(7ra/2) is exactly as expected from the analytic continuation: 

<->■ (—loj + 0^)““ = w““[cos(7ra/2) + isin(7rQ!/2)]. (A4) 

We can also recover the real part of by making use of the Kramers-Kronig transform (an application of the spectral 
representation), 


Rey^(a; + i0+) 


Jq tt a;'2 - 


(AS) 


where we have again made use of the fact that p{uj) is odd. Vj denotes the principal value of the integral. Using an 
infrared cutoff ujq, we find that the real part corresponding to a spectral power law tail l/u;“ is 


r 

Rey^(a; + i0+) 

J ijj 


dw' 2 cli' 1 


U)0 

cot(7rQ:/2) 


TT 


Wn 


,2n 


n—1 


(A6) 

(A7) 


where cot x = cos(a;)/ sin(a:), and bn = l/(a — 2n). The first term is independent of wg and precisely yields the answer 
expected from Eq. (A4), whereas the remaining series again decouples because it only contains even powers of l/w. 

In the above, we have assumed that p{uj) decays to zero as a; —>■ oo. If this is not the case, one can analytically 
continue the terms with positive powers of w using Eq. (A4), and apply the above procedure to the modified sus¬ 
ceptibility, y — y — i5y°°, which vanishes as w —)■ oo. We refer the reader to the main text for a detailed discussion 
regarding the subtractions y, Sx°°- We have thus shown that the asymptotic expansions of the retarded and Euclidean 
susceptibilities are precisely related by a termwise analytic continuation. 

Our results can be explicitly checked for the infrared hxed point of the 0{N) model in the N —>■ oo limit, as is 
discussed in the main text and in the next appendix. In the context of the charge conductivity in 2-l-lD, this was 
done for the 0{N —)■ oo) model and the Dirac GET in Ref. 6, and for a wide class of conformal QC theories without 
quasiparticles using AdS/CFT in Refs. 6, 14, and 17. 


1. Even powers in the spectral function 


We now turn to the more subtle case of even powers in p{u}). From the basic formula for analytic continuation, 
Eq. {A4), we see that even powers in imaginary frequencies contribute only to the real part of y^. E.g. 1/U^ —>■ —l/w^, 
which does not contribute to p = Imy^. To understand what asymptotic Euclidean term can give rise to a l/uP 
scaling for p, we turn to the spectral representation, Eq. (A2), 


dcj' 2uj' 11 / ^ 









TT ^ 
n—2 


(A8) 


(A9) 
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where ujq is the infrared cutoff we have employed above, and r„ = (—l)"/(n — 1). The first term decouples from the 
series due to the presence of the logarithm; also note that there is no term in the series. We thus see that a 
logarithm in y® gives rise to a l/oj^ term in p. This is fully consistent with the direct analytic continuation of the 
term with the logarithm: 


Inf^) + In (4) 

\^o) tn-»i^+to+ _ V‘^0/ 

TTU)'^ 



(AlO) 


for w > 0. Analogous results hold for higher powers of l/w, e.g. l/w^ obtains from ln(ff^/a;Q), etc. We note 

that such logarithms do not appear in the leading expansion of the response functions of the Wilson-Fisher QC points 
studied in this paper. 


Appendix B: Scalar susceptibility of the 0{N) model at large N 


We give details regarding the asymptotics and sum rule for the scalar susceptibility Xs{k) = {Og{—k)Og{k))T of the 
0(A) Wilson-Fisher fixed point in the large N limit. Using the 0(A) non-linear sigma model defined in the main 
text, it can be shown^ that for A —>■ oo the susceptibility is given by 


Xs{k) 


A 

nM’ 


(Bl) 


with the scalar polarization function 


nE(i,,„,k) = T^ 

Vji 


d'^q 

(27r)‘^ 


X® {iujn + k -I- q)xn {ivn, q) • 


(B2) 


The sum/integral involves 


xlbik) = 


the order parameter susceptibility (propagator). At zero external spatial momentum, we find 

d'^k l-|-2nB(efc) 1 


n^(za;„, 0 ) = J 


44 +^; 


= K, 


{2t:Y efe 
, Tde (e^ - ^ 

JmT 4e^ -I- a;„ 


(B3) 


(B4) 

(B5) 


where e\ = k^ -I- m|., ubY) = l/[exp(e/T) — 1] is the Bose-Einstein distribution, Kd = j 2 ^ = 

(T{dl2) is the normalized area of the unit sphere Sd-i, and r( 2 ) is the Gamma function. 


1. Asymptotics 


We find the following uJn^ T expansion for 11®: 


n®(ia;„,0) = 1 + 

\ 


{d + l)T{d/2)r{-{d + l)/2)m: 


d-1 

T 


40r 


poo 

/ de(e^ - my)^'^“^)/^2nB(e) 
J 


UJ^ 


rpd-\-l 

. ,d+l 


+ ■ 


, (B6) 


















FIG. 2. Thermal mass rriT/T = Qd, as well as expansion coefficients ao,ag appearing in Xs(l^l T) as a function of the 
spatial dimension d. Interestingly, ag vanishes in 2 + ID. 


where the dimensionless coefficients are given by: 

Kdn sec{dn/2) 

“0 = ; 

KdQd ~ 2)’’" sec(d7r/2) 


Go 


2d-l 


QjT = 


I r(d/2)r(-(d + i)/2) 

Oo 


/ OO 

de(£2 _ l)(rf-2)/2^2 

nsirriTe) 


(B7) 

(B8) 

(B9) 


0d, aoj Og are plotted in Fig. 2 for I < d < 3. The coefficient of the term arising from the stress tensor, ot, is positive 
for 1 < d < 3. We hnd that the coefficient of in Eq. (B6) vanishes exactly once the value of the thermal mass 

rriT is used (see Sec.B2). This is as expected since the QC point does not have a 0(Ai)-invariant scalar with scaling 
dimension d— I, and thus a term in Eq. (B6) would be at odds with the OgOg OPE. An analogous cancellation 

was found® to occur for the conductivity of the 0{N) model. Eq. (B6) then leads to an asymptotic expansion of xf: 


X^{k) 



+ Gt 


'Y'D 

w 



(BIO) 


where we have reinstated the the full dependence on j/cp = o;^ + k^. In two spatial dimensions the above simplifies 
to: 


X^{k) 


-8Al|A:|® 



28C(3) T® 

Stt |/c|® 


(Bll) 


Interestingly, from Eq. (B8) we note that Gg vanishes exactly in d = 2 spatial dimensions (but not when d yf 2). This 
stems from the fact that the OPE coefficient of Og in the OgOg OPE, Cggg, vanishes exactly in d = 2 and N = oo. 
This was previously noted in Ref. 23, where the expansion Eq. (BIO) was also given. As mentioned in the main body, 
we do not expect that Cggg vanishes at finite N. In agreement with this expectation, the OPE coefficient in the Ising 
case {N = I) was recently computed^^ by means of Monte Carlo simulations and found to be finite. It would be 
interesting to compare this new result with conformal bootstrap or with a 1/A^ expansion on the field theory side. 


2. Thermal mass 


Interestingly, by imposing the vanishing of forbidden terms in the asymptotics of a 2-point function such as Xs, we 
can determine the value of the thermal mass, mr = QdT- Indeed, reverting back to the integration variable k in 
Eq. (B6), we find that setting [• • • ] =0 leads to the following integral equation for 0^: 


d^k UBjek) 
{2nY Cfc 


QdT'^-^Xd+i + 


(BI2) 
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where + {Qd.T)^, and we introduced the dimensionless constant Xd+i = 2r((3 — d)/2)/[(47r)(‘^+^^/^(d— 1)]. 

This agrees exactly with the equation obtained by requiring that ipa^Pa = 1 in the non-linear sigma model.In D = 3, 
this equation can be exactly solved:® 0^ = 21n(/5, where <p = [1 + \fh')l2 is the golden ratio. 


3. Response function and sum rule 


The scalar response function is defined as follows: 


7?.s(w,k) 


zw - 0+ ’ 


(B13) 


where Xs is the susceptibility. We shall set the spatial momentum k to zero for simplicity. Our goal is to explicitly 
verify the sum rule for TLs- Using the general result given in the main text, Eq. (6), we have at iV = oo: 

dw Re [n,{uj) - 7^f=°(w)] = 0, (B14) 

which is precisely Eq. (10). Using the asymptotics obtained above, we find that = {N/ao){u}/i)^~^, which 

becomes frequency independent for D = 3. After subtracting this leading asymptotic behavior, the integrand decays 
sufficiently fast for the sum rule to be well-defined. Indeed, using Eq. (BIO), we have for oj T: 



7^«(a;)-7^f=°(^c) = 


iagN/ao)T^ 


(B15) 


which comes from the relevant scalar Og, with scaling dimension A^~°° = 2 for 2 < U < 4. 

We now provide some details regarding the numerical verification of the sum rule. In general, the real part of the 
response is given by: 


ReTZsioj) = -7rxf{0)S{uj) +-Imxfioj) ■ 

LJ 

At N = oo, we have®’® xfi^) = so that the spectral function for xf becomes 

^Imn«(u;) 


Imxf(w) = N 


in^Hr 


(B16) 


(B17) 


where |n^(a;)| is the complex norm. The imaginary and real parts of n'^(tLi) at w > 0 can be obtained from Eq. (B5): 

Imn^(w) = 0(a; - 2 tot)^ - ^t] [1 + 2nB(w/2)]; (B18) 

Aoj 

Ren^ {uj) = KdV [ de(e2-m|)(‘^-2)/2l^^!^, (B19) 

where 0(a:) is the step function, and U/ denotes the integral’s principal value. We note that Imxf'(a;) and Re77.s(w) 
vanish for 0 < a; < 2mT, where tot is the “thermal mass” of the quasiparticles at = oo, Eq. (B3). This hard 
gap behavior is an artefact of the N = oo limit. We further note that Re 11^ has a logarithmic divergence at the 
threshold oj = 2mT- The resulting numerical plots for Re77s(w) are given in Fig. 3a. Ref. 40 has computed Imxs(w, k) 
at frequency and momentum k, where additional subtleties in the integrals emerge. The calculations of Ref. 40 can 
be used to verify the momentum-dependent sum rules, a task we leave for future investigation. We finally note that 
Re77s has a delta function at a; = 0, Eq. (B16), the weight of which is given by 


-7rxf(0) 


Ntt 

n«(o) ’ 


(B20) 


which is plotted in Fig. 3b as a function of D. For 77 = 3, we obtain N x 53.04194. 

The integral for the sum rule Eq. (B14) can be divided into three parts: the w = 0 delta function, the interval 
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FIG. 3. a) Scalar response of the 0{N) model = oo) in various spacetime dimensions D. The dashed vertical arrow 
represents a delta function S{u}). b) Weight of this delta function versus D. 


(0, 2tot) where Re7?.s vanishes, and (2mT,+oo): 

|•2mT 

I = - 2 Xf (0) + dc. Re [O - <=°(u;)] 

For D = 3, we numerically find: 


poo 

/ dwRe[7^s(a;)-7^J=°(w)] . 

J 2mT 


(B21) 


- 2 Xf(0) =+26.520972; 


p2mT 

- / da;Re7^f=°(w) = -15.398778; 

Jo 

pOO 

/ dwRe [77,(a;)-77f=°(a;)] =-11.122193. 

J 2rm- 


(B22) 


The terms sum to X = 1.9 x 10“^. We can attribute the deviation from zero to the numerical uncertainty in our 
evaluation of the integrals. (We emphasize that care must be used to treat the logarithmic vanishing of ReTT-s at the 
threshold ui = 2mr.) The fact that three numbers of order 10^ cancel to within 10“^ provides an excellent check 
of the sum rule. We have also verified the sum rule for other spacetime dimensions D, and have found that it is 
respected, in accordance with the proof given in the main body. 
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